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Abstract 

We solve exactly the classical non-relativistic Landau-Lifshitz equa- 
tions of motion for a charged particle moving in a Coulomb poten- 
tial, including radiation damping. The general solution involves the 
Painleve transcendent of type IL It confirms our physical intuition 
that a negatively charged classical particle will spiral into the nucleus, 
supporting the the validity of the Landau-Lifshitz equation. 
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1 Introduction 



A corner stone of theoretical physics is the exact solution of the motion 
of a particle moving in an inverse square law force. The orbits are conic 
sections: ellipses or hyperbolae depending on initial conditions. The original 
application was to the motion of planets around the SunpLj . Later the same 
problem was found to arise in the Rutherford scattering of alpha particles 
and in the classical model of the atom. 

The orbits of charged particles in a Coulomb potential cannot be conies 
exactly, as it is a fundamental principle of electrodynamics that all charged 
particles must radiate when accelerated [21 E]. The radiation carries away 
energy, and therefore acts as a dissipative force, changing the equation of 
motion. Can we still find an exact solution for the motion of the particle in 
a Coulomb field, after taking account of radiation reaction? 

Deriving the correct equation of motion for a charged particle, including 
this radiation damping, is not a simple matter. The problem is that radiation 
reaction is the force exerted on the particle by its own electromagnetic field; 
a straightforward application of the Lorentz force law will give an infinite 
force in the case of a point particle. Dirac [1] found a way through this 
minefield of divergences and deduced an equation of motion including the 
radiation reaction. A key point was that the divergences are removed by a 
renormalization of the mass of the charged particle. 

This Lorentz-Dirac equation of motion is a third order ODE, as the radi- 
ation reaction force is proprortional to the derivative of acceleration. Typical 
initial conditions will give unphysical solutions that 'runaway': the energy 
grows without bound instead of decaying. Thus Dirac's work, although a 
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major step forward, cannot be the final word on the equation of motion of 
charged particles. 

Spohn [3, 5j showed that these unphysical runaway solutions can be 
avoided if the initial data lie in a 'critical manifold'; i.e., if the initial condi- 
tions are fine-tuned to avoid the runaway unphysical solutions. This turns 
out to be the same as treating the force due to the radiation as a first order 
correction. We get this way a second order equation with physically sensi- 
ble solutions. Although without the modern understanding, this equation 
of motion for a radiating particle were given first in the classic text of Lan- 
dau and Lifshitz^. Therefore these are known now as Landau-Lifshitz (LL) 
equations of motion. See Ref. for a physical argument in support of the 
LL equations. 

There is no unanimity yet that these are the exact equations of motion of 
a radiating charged particle [8] . In addition to experimental tests, we need to 
verify their theoretical consistency. As an example, it should not be possible 
for a negatively charged particle to orbit a nucleus in an elliptical orbit: it 
should plunge into the nucleus as the radiation it emits carries away energy 
and angular momentum. It is important to verify the physical correctness of 
the LL equations by checking that its solutions have this property. 

In this paper, we solve exactly the non-relativistic Landau-Lifshitz equa- 
tions in a Coulomb field ; the general solution is in terms of [9l [10] Painleve 
transcendents of type II. The same differential equation , with different initial 
conditions, appears in several other physical problems, such as the Tracy- 
Widom law for random matrix eigenvalues [TT| [12] . Our solution turns out 
to have the correct asymptotic properties: the orbit of a negatively charged 
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particle does spiral in towards the nucleus. 

The Lorentz-Dirac equation of motion has been studied in the Coulomb 
field. It has a complicated, unphysical behavior, and no general exact so- 
lution is known. For example, in an attractive Coulomb potential there are 
solutions that accelerate away to infinity. See Ref. |3j Section 6-15, [13] The 
approximate, numerical treatment of the radiative reaction in Ref. [H] is 
closer to our physical results . 

Due to quantum effects, the LL equations cannot be the right description 
at the short distances characteristic of atoms. Our solution might still be 
an approximate description of an electron with a large principal quantum 
number in an atom. It should also describe an alpha particle scattered by a 
nucleus and an electron (or positron) emitted by a nucleus, all of whose mo- 
tion is affected by the radiation emitted. Relativistic corrections will become 
important as the velocity of the particle grows; we are currently investigating 
the exact solvability of the relativistic LL equations in a Coulomb field. 

The LL equations have already been solved in a constant electric field 
and in a constant magnetic field[T5]. We hope that more physically realistic 
situations will open up to study using the techniques we describe here. In 
another paper [TB] we have proposed a canonical formulation and a quantum 
wave equation for dissipative systems of a particular type. The case we study 
here happens to be of this type, so we hope that a quantum treatment of a 
radiating electron in an atom along these lines is also possible. This might 
allow us to go beyond perturbation theory in the calculation of line-widths 
of a hydrogenic atom. 

The analogous problem in General Relativity of a star being captured 
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by a blackhole, its energy and angular momentum being carried away by 
gravitational radiation is of great importance in connection with the LIGO 
project to detect gravitational waves. We hope that our solution of the much 
simpler electro dynamic problem will help in understanding this well. 



2 The Landau-Lifshitz Equations 



The LL equation of motion of a radiating charged particle in an electrostatic 
field is [6], 
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q and m being the charge and mass of the particle respectively. The dissipa- 
tion parameter r has units of time; for the electron it would be the (2/3 ) of 
the classical electron radius divided by c. 

In the non-relativistic limit, it is much simpler: 

llv + rVtT + Vf/ = 0, v = | (4) 

where U is ^ times the electrostatic potential. For a central potential, 

VU = rUr, Ur = ^ (5) 
ar 

^ [v + rrUr] + rUr = 0. (6) 



Taking the cross product with the position vector gives, with L = r x v, 

^ = -r^L. (7) 

Thus the direction of angular momentum is preserved. If the initial con- 
ditions arc such that L 7^ 0, the motion will lie in the plane normal to this 
vector. If L = initially, it remains zero and the motion is along a straight- 
line passing through the center of the potential. 

Using the standard identities 

r2 



^v^, + - (8) 



dr ^ d ^ \ , ^ d ^ ^ 

— ^ r.v, — r = - [v - VrV] , r.— r = (9) 

dt dt r dt 

dvr L"^ ^ dy , , 

we get the system of ODE 

dr d , dL Ur ^ 

3 The Coulomb Potential 

For the Coulomb potential U — ^,Ur — and 

dL krL 



dt r3 

Thus the centrifugal force is a total time derivative: 

1 d 



This coincidence allows to write the radial force equation as 



~dt 



r 1 j2 
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Put 



to write this as 



kr 



1 



Z = Vr 



dL kr dz k dr 



kr 



dt ' dt 



dt 2kT 



(15) 



(16) 



Using the fact that this is an autonomous system ( i.e., t does not appear 
exphcitly) we can ehminate dt, to get a system of two ODEs, 



dL T ^ dr 1 Ir, 

— — T L -\- T Z -\- T 

dz r ^ dz 2A;2r^ k 



(17) 



We note as an aside that in the case of purely radial motion, L = this 
reduces to a Riccati equation for p = -: 



dp 

dz 



k 



rp 



(18) 



This can solved in terms of Airy functions. 

Returning to the general case, we can rewrite the above system as a single 
second order ODE: 

(19) 



dz^ 



2P k 



Up to scaling, this is the particular case with a = of the Painleve II equation 
( See Ref. 0, page 345 ) 

d^u 



dx^ 



2u + xu + a 



Defining constants 



we have the solution 



, a 
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L = au(bz). 



(22) 



When A; < 0, as for an attractive Coulomb potential, u is purely imaginary 
and the independent variable x = bz is real, u also depends on a complex 
'modular' parameter s that is determined by the initial conditions (TO]. 
The identities 

odO dz _ k 

^-'Tt' ^^^^ 

allow us to determine the polar angle: 

By a quadrature and a differentiation of the Painleve transcendent, r and 9 
are both found as functions of the parameter z, determining the orbit : 

r{z) = rg, eiz) =Oi + l r u{bz)dz. (25) 

4 Examples 

4.1 A Decaying Orbit 

To plot orbits, another form of the equations is sometimes more convenient. 
Define y = L?' and change to 9 as the independent variable: 

dL kr dz L z kr 



d9 r ' d9 2kT L Lr 
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(26) 



dL' d 
2kr 



d9 r' d9 

to get the third order ODE: 



L 
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2kT 
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d^y dy 2A;^r „ 
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The orbit is then given by 



1 



1 dy/y 



(29) 



r{9) kr d9 

We can find the orbit by numerically integrating the above third order ODE 
This was done before it was realized that the equation can be solved analyt- 
ically. It is still useful as a way to plot a slowly decaying orbit: 
y 




4.2 A Capture Orbit 

We want a solutioij^l for which r = ^ > . Note that ^ < 0. If r(^) 

vanishes it must be a simple zero, since near a zero ^ ~ 1- Thus r ~ z — zq, 

L C{z — zq) where C and zq are real constants of integration. 

If a particle approaches from infinity with angular momentum Li and 

^ For simplicity, we use in this section the natural units of the problem, with |fc| = r = 1. 
Normal units can be restored by dimensional analysis. 
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radial velocity Vi, 



Since ^ ~ tiir^, we get 



zi = vi + -L^. 



(30) 



(31) 



Vi[Z - Zi) 

Thus zi is a simple zero of ^/'(z). 

In other words, to describe the capture of an incoming charged particle 
by an attractive Coulomb potential, we just have to solve the Painleve II 
equation with the initial conditions 



zi = vi + ^Ll, u{zi) = ^=Li, u\zi) = 



(32) 

and evolve to a point zq < z\ at which ^(^^o) = 0. The complete orbit 
corresponds to the finite range zq < z < zi in the parameter z. We plot an 
example, obtained by numerical calculations, of such a capture orbit below. 




-3.0 -2.5 -2.0 -1.5 -1.0 -0.5 



The velocity blows up as ^ for small r. From a finite distance the particle 
is captured in a finite time. 
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Other cases can be worked out similarly. The formulation of the Painleve 
equation in terms of the Riemann-Hilbert problem and Isomonodromy [T2] 
give powerful techniques to study our solution. In particular, there are a pair 
of conserved quantities that replace energy and angular momentum in this 
dissipative but integrable system. Also, we can generalize the Rutherford 
formula for Coulomb scattering to include radiation. We hope to return to 
such a detailed analysis in a longer paper. 
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